KUNS-1948 
KANAZAWA-04-19 



Democratic mass matrices induced by strong gauge 
dynamics and large mixing angles for leptons 

Tatsuo Kobayashi"^, Hiroaki Shirano^^ and Haruhiko Terao^^ 

"Departmant of Physics, Kyoto University, 
Kyoto 606, Japan 

''Institute for Theoretical Physics, Kanazawa University, 
Kanazawa 920-1192, Japan 

Abstract 

We consider dynamical realization of the democratic type Yukawa cou- 
pling matrices as the Pendelton-Ross infrared fixed points. Such fixed 
points of the Yukawa couplings become possible by introducing many 
Higgs fields, which are made superheavy but one massless mode. Ex- 
plicitly, we consider a strongly coupled GUT based on SU{5) x SU{5), 
where rapid convergence to the infrared fixed point generates sufficiently 
large mass hierarchy for quarks and leptons. Especially, it is found that 
the remarkable difference between mixing angles in the quark and lepton 
sectors may be explained as a simple dynamical consequence. We also 
discuss a possible scenario leading to the realistic mass spectra and mix- 
ing angles for quarks and leptons. In this scheme, the Yukawa couplings 
not only for top but also for bottom appear close to their quasi-fixed 
points at low energy and, therefore, tan /? should be large. 
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1 Introduction 



Yukawa couplings are not direct observables, since they change depending on the choice of 
flavor basis. What we know in principle are the hierarchical masses of quarks and leptons, 
and the mixing matrices. The recent neutrino oscillation experiments [1, 2, 3, 4] indicate 
strongly that neutrinos are massive and lepton flavors are mixed. It is remarkable that the 
mixing angles (sin^ 2^sun ~ 0.84 and sin^ 2^atm ~ 1-0) are as large as nearly bi-maximal, 
which are in sharp contrast with the small mixing angles among quarks. 

There have been proposed various phenomenological models explaining origin of the 
Yukawa couplings accommodating these features by assuming certain flavor symmetries. 
[5]. The most popular one would be the Froggatt-Nielsen mechanism [6] based on abelian 
flavor symmetries. In this line of thought, hierarchical structures among the Yukawa 
couplings are nicely explained. However the abelian symmetry itself cannot flx 0(1) 
ambiguity of parameters. Therefore structure of the mixing angles in the lepton sector 
as well as the neutrino masses are explained no more than that those parameters are not 
hierarchical. Also the flavor differences are just input in the abelian flavor charges apriori. 

Contrary to this, non-abelian flavor symmetries can impose flrm relations among some 
elements of the Yukawa matrices. Therefore various non-abelian symmetries have been 
considered to explain the mixing angles in the lepton sector as well as the hierarchical 
masses simultaneously [7]. The democratic Yukawa matrices have been an attractive 
possibility of this kind [8, 9, 10, 11, 12]. The rigid democratic matrix, in which all 
elements are identical , may be explained by assuming e.g. S3 flavor symmetries for 
three generations of quarks and leptons. It is amusing that one of the mass eigenvalues 
appears much larger than others, though there is no apriori difference among flavors. 
The 5*3 flavor symmetry is supposed to be broken shghtly to give rise to small masses 
for the flrst and second generations. It has been known for some time that Fritzsch 
type mass matrices for quarks are realized after introducing a simple form of 6*3 breaking 
parameters [8]. Moreover the lepton sector mixing angles uncovered by neutrino oscillation 
experiments may be explained in the democratic framework rather well [9, 10]. 

However no comprehensive explanation has been given to the origin of the flavor 
symmetry breaking terms in this framework. Moreover the 5*3 flavor symmetry allows the 
mass matrix proportional to an identity matrix other than the democratic type matrix 
for the neutrino sector. The mixing matrix for the lepton sector can be obtained, only 
if we assume the neutrino mass matrix to be almost diagonal. Namely the democratic 
part must be small with some other reasons than the flavor symmetry [9, 10]. The same 
problem remains in the see-saw context, where both of Yukawa couplings and mass matrix 
of the right-handed neutrinos should be almost diagonal. 

On the other hand, we face with a more stringent problem in considering the grand 
unifled theories. In e.g. the SU(5) GUT, the matters belong to a 10 and a 5. However 
S3 symmetry allows an identity matrix for the Yukawa couplings among lO's, which must 
be strictly forbidden for the quark mass hierarchy [11]. Thus the phenomenologically 

^Phenomenologically the so-called extended democratic mass matrices [13] are equally considerable. 
However we restrict ourselves to the rigid cases in this paper. 
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postulated forms for the democratic Yukawa matrices are not explained only by flavor 

symmetries. 

Indeed it is not a unique way to introduce some flavor symmetries and their small 
breakings in order to explain hierarchy among various couplings. For example, large mass 
hierarchy may be generated by overlapping of the wave functions in extra dimensions. 
Renormalization with large anomalous dimensions induced by strong dynamics also can 
generate mass hierarchy [14]. Similarly power-law running of the Yukawa couplings in the 
extra dimensions has been also utilized [15]. 

In this paper we consider the models in which the democratic type of Yukawa matrices 
are reahzed as infrared attractive fixed point couplings [16, 17, 18, 15]. It has been known 
for some time that ratio of the Yukawa couplings and the gauge coupling approach rapidly 
to the so-called Pendelton-Ross fixed points (PRFP) towards infrared [16]. If all elements 
of the Yukawa coupling matrix have an identical PRFP, then the democratic type matrix 
may be achieved dynamically. In order to make this possible, we need to introduce many 
Higgs fields. Also a special kind of mass terms for the multi Higgs fields are assumed so 
that they obtain masses of the GUT scale except for one massless mode, which is identified 
with Higgs in the low energy theory. 

Indeed this idea is not new. Abel and King[18] have already considered several years 
ago except for the neutrino sector. In this paper, however, we stress that this mechanism 
is free from the difficulty in realizing democratic Yukawa matrices by the fiavor symmetry. 
Moreover we are going to show that the typical difference in the mixing angles between 
quark and lepton sectors may be explained as a simple dynamical consequence in this 
framework. Our mechanism predicts very small neutrino Yukawa couplings, and therefore 
relatively light right handed neutrinos in the see-saw models. 

Explicitly, we treat a super symmetric GUT model based on SU{5) x SU{5) gauge 
group, in which the doublet-triplet problem is solved nicely [19, 20]. There one SU{5) 
gauge interaction coupled with quarks and leptons is assumed to be strong and another is 
weak just as in the conventional SU(5) GUT. Then renormalization group (RG) from the 
Planck scale to the GUT scale shows that the Yukawa couplings in the GUT model are 
aligned in a very good accuracy due to strong dynamics. After diagonalizing the Yukawa 
matrix, only one coupling is found to be much larger than others. Further we will show 
that the realistic Yukawa coupling matrices may be obtained within a simple setup. 

In this kind of scenarios, all of the Yukawa couplings, except for the neutrino sector, 
are given to be quite large at high energy. Because the Yukawa coupling and the gauge 
coupling are of the same order at the PRFP. However heavy top quark mass (178GeV) 
indicates that top Yukawa couphng is close to the so-called quasi- fixed point [21]. This 
means that top Yukawa coupling can be fairly large at the GUT scale, which is a favorable 
point of our scenario. However Yukawa couplings not only of top quark but also of bottom 
quark and tau lepton are all large. Therefore large tan/3 is predicted in this kind of 
scenarios. 

The article is organized as follows. In section 2 we give a brief summary of the 
democratic type of mass matrices and the 5*3 flavor symmetry. In section 3 we present 
the general ideas leading to the democratic Yukawa matrices at low energy by strongly 
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coupled gauge dynamics. There also the superpotential for the multi Higgs fields and 
their mass spectra are discussed. In section 4 we give an explicit model based on the 
SU{5) X SU{5) GUT. Then it is seen that sufficiently large hierarchy in couplings can be 
generated by RG from the Planck scale to the GUT scale. We also discuss predictions 
of this kind of scenarios and consistency with particle masses. There we work out the 
Yukawa couplings for top, bottom and tau obtained at low energy scale in the explicit 
model. The neutrino Yukawa couplings are discussed in section 5. There it is shown that 
the large mixing angles among leptons may be realized as a dynamical consequence, if 
the right-handed neutrinos have new Yukawa couplings to some other fields. In section 6 
we consider minute structure of the democratic mass matrices producing reahstic masses 
and mixing angles. Finally we devote section 7 to conclusions and discussions including 
comments on the extra dimensional setup and also on the soft supersymmetry breaking 
parameters. 



2 Democratic mass matrices and 6*3 flavor symmetry 

In this section we give a brief review of the democratic mass matrices for quarks and 
leptons from the view point of 5*3 fiavor symmetry. The democratic mass matrices give 
a phenomenologically successful description. However it is found to be necessary to con- 
strain the parameters allowed by the symmetry further. Namely the symmetry is not 
sufficient to give viable mass matrices. Indeed some attempts have been done in order to 
give comprehensive grounds to the viable democratic mass matrices so far [12]. 
The rigid democratic matrix J is diagonalized as 
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by the diagonalization matrix 



A 



1/V2 
-1/V2 




-2/^/6 



(2) 



For the democratic mass matrix, two vanishing eigenvalues are regarded as masses for the 
first and the second generations in the first approximation. Their masses are given by 
small deviations from the rigid matrix. It would be said that the interesting feature of 
this approach is the idea of flavor democracy, namely that there is no difference between 
flavors apriori. 

In the standard model (SM), the democratic mass matrices are realized by assuming 
fiavor symmetry in each kind of the SM matter fields. Hereafter we consider the minimal 
supersymmetric case, the MSSM. The superpotential of the MSSM is given as 



R 



(3) 
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where i,j — (1, 2, 3) represent generations. Each of Qi,Ui,di, Li, is assigned to a three 
dimensional representation of a distinct S3. Then the Lagrangian is invariant under 
permutation of any set of the matter fields. 

Let us start with mass matrices for quarks. It has been known that the realistic mass 
matrices are obtained by introducing small breaking parameters eu{d) ^ Su{d) -C 1 of the 
flavor symmetry as follows [8]. 



J + 





-eq 













eq 





I 











= A 



'l/3e, 
(2/3)5, 



2/3e, ^ 



(4) 



[-^2/3eq -{V2/3)6q 1 + {1/3)6 J 



where q = u, d. The matrix in the 2nd line is known as the Fritzsch type, which offers us 
a phenomenologically good representation. 

The mass eigenvalues of these matrices are found to be 



m? : m% : m% ~ : 5„ : 1. 



(5) 



Therefore ratio of quark masses can be represented by assuming the breaking parameters 
of the following order; 

lo-^ 5.-10-2, 

-2 x _ , y^) 



10" 



The absolute values of the top and bottom quark masses are not determined by the flavor 
symmetry and, therefore, must be parameterized phenomenologically. In our dynamical 
mechanism, wc can predict these absolute values by evaluating the Yukawa couplings at 
low energy. This will be discussed later. 

The mixing angles turn out to be small, since the diagonalization matrix for and 
Md are both close to A. Quantitatively the CKM matrix is given in a good approximation 

by 

/ 1 OieJSa) Oiea) 



UcKM 



\ 



1 



(7) 



0{eJ5a) 1 
0(e,) 0{5a) 

It is found that the above choice of parameters may reproduce the realistic mixing angles 
also. 

Next let us go into the lepton sector. The mass matrix for the charged leptons are 
given in a similar fashion; 



Me OC J + 









(8) 
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Here, if we take the flavor symmetry breaking parameters as 

ee - 10-2, 5, - (9) 

then ratio of the charged lepton masses are parameterized well. Again the absolute value 
is undetermined in this approach, though it will be determined in our scheme. 

However situation for the neutrino mass matrix is somewhat different. The fla- 
vor symmetry allow the couplings nij given by Eq. (3) proportional to an identity also. 
Therefore the mass matrix is given as 

/ \ 

oc / + r J + ey , (10) 
V 5, I 

where / denotes an identity matrix. 

It is noted that this mass matrix may explain the mixing angles in the lepton sector 
nicely, if the parameter r in Eq. (10) is much smaller than 1. For vanishing r, the lepton 
mixing matrix (or the MNS matrix) is just given by A'- (matrix A is given explicitly in 
Eq. (2)). Therefore we obtain the three mixing angles as 

sin^ 2^12 ~ 1, sin^ 2^23 ~ 0.94, [7^3 = 0. (11) 

This should be compared with observation by the recent neutrino experiments [1, 2, 3, 4], 
which tells us 

sin^ 2^12 ~ 0.84, sin^ 2^23 ~ 1-0, If/esI < 0.23. (12) 

Thus the lepton mixing angles are found to be almost explained already without parameter 
r in Eq. (10). In other words, a nearly diagonal neutrino mass matrix is favorable including 
the flavor symmetry breaking part from phenomenological point of view. However we do 
not flnd any reasons either why r is so small, or why the flavor symmetry breaking 
parameters are input only at diagonal elements. 

This fine-tuning problem remains or becomes more curious in the see-saw mechanism. 
The supcrpotential for the neutrino sector is given by 

W^Y^^.LiV^H- + \mr,^v,v^, (13) 

where z/j denote the right-handed neutrinos. Here we need to assume that these right- 
handed neutrinos belong to the same representation of 5*3 as the lepton doublets Lj. If we 
assume a distinct group for the right-handed neutrinos like other fields, the neutrino 
Yukawa coupling matrix Y" is restricted to the democratic form. Then the mixing angles 
cannot be large. The neutrino Yukawa matrix and the right-handed neutrino mass matrix 
are now parameterized as 

= yl{I + TJ) + Ay^ (14) 
Mr = MRo{I + r'J) + AMr. (15) 
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In should be noted that both of V and Mr should be nearly diagonal in order to have 
large mixing angles . Therefore the free parameters r and r' are constrained to be small, 
which is not explained by the flavor symmetry. 

Moreover this sort of flne-tuning turns out to be required much stronger in the GUT 
models. In the SU(5) GUT, the quark and lepton fields in one generation are combined 
into a 10 + 5 + 1 representation. The superpotential of their Yukawa interactions is given 

by 

W = Yi; 10a0,i/(5) + 10,5%i/(5*) + Y^"- 5*a,if(5). (16) 

We note that Yukawa couplings satisfy Y'^ = {Y'^)'^ at the GUT scale. Then symmetry 
for lOj allows the Yukawa coupling matrix of F " = y^iJ + r"I) in general. However the 
mass hierarchy of up-sector quarks cannot be obtained unless r" is suppressed to 0(10^^). 

Thus the approach of democratic mass matrix is an attractive possibility phenomeno- 
logically, but the fiavor symmetries do not support origin of the assumed Yukawa cou- 
plings. In the next section, we consider realization of the democratic Yukawa matrices in 
a dynamical fashion without recourse to any flavor symmetries. There it will be seen that 
we are not troubled with the flne-tuning problems any more. 



3 Flavor democracy by strong gauge dynamics 

Our basic assumption is that the Yukawa couplings are not hierarchical nor aligned, but 
may be even somewhat anarchy at the fundamental scale. The purpose of this section 
is to show that the Yukawa couplings can be aligned to the rigid democratic form quite 
minutely at lower energy scale owing to strong gauge dynamics. The realization of realistic 
matrices are discussed in the later section. Also we consider such a mechanism in GUT 
models, since the problem of fiavor symmetry is very severe especially in GUT. 



3.1 IR fixed point for a single Yukawa coupling 

Before going into the democratic Yukawa matrix, let us consider the IR fixed point a la 
Pendelton and Ross [16] in the cases with a single Yukawa coupling. In general, the RG 
equations for the gauge coupling g and the Yukawa coupling y are given at one-loop level 

as 

/x-j-^ — {aay — cag)ay^ (18) 
(jjfjj 

where we defined ag = g^/Sn"^ and ay = ||/p/87r^. The coefficients depend on the fields 
contents and the Yukawa interaction, but note that a and c are positive constants. Then 
the ratio x — ay/ ag satisfies the RG equation, 

ji— — [ax — (c — h)\ agX. (19) 
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11-^ = (c - b)agAx. (20) 



This beta function has a non-trivial fixed point at x — x* — {c — b)/a, which is IR 
attractive. 

The convergence behavior of the RG flows around the flxcd point may be seen by hnear 
analysis. The deviation from the fixed point, Ax = x — x*, is subject to the equation 
given by 

dAx 
d/i 

If the gauge coupling does not change rapidly, then the fiows of x around the IR fixed 
point may be evaluated as 

Axin) ^ (^j Aa;(A). (21) 

We see that the conditions for strong convergence as follows; (i) The gauge couphng is 
large, (ii) The coefficient (c^b) is positive. Also larger (c — b) makes convergence stronger. 
The second condition means that the coefficient b should not be very large, therefore, that 
asymptotically free gauge theories with rapidly running gauge couplings arc excluded. ^ 
Such a fixed point appears in gauge theories beyond four space-time dimensions as 
well[15]. In the dimensions of 4 -|- 5, the RG equations for the gauge couphng and the 
Yukawa coupling at one-loop are given as 

//^ = 5ag-bag^ (22) 

/x-y-^ = Say + {aay — ca^)ay, (23) 

CLlJj 

where ag and ay denote dimensionless couplings. Then the RG equation for x — oty/oig 
is found to be identical to Eq. (19). We should note some particular features in the 
extra-dimensional theories. First the running law of these couplings is not logarithmic 
but power with respect to the renormalization scale. Also the gauge coupling becomes 
strong at ultraviolet irrespectively of the sign of b. With these reasons, convergence to 
the IR fixed point is found to be very strong in general. 



3.2 Models with 3 flavors 

If we extend the Yukawa interaction to multi-flavors naively, then the IR fixed points 
for the Yukawa matrices may exist but are proportional to an identity matrix. Thus 
the democratic Yukawa coupling matrix cannot be obtained by renormalization effect in 
MSSM. However it turns out to be possible once we admit multi Higgs fields. 

First let us consider only the up-type Yukawa couplings in the superpotential given 
by Eq. (3) or by Eq. (16). We introduce 9 elementary higgs superfields Hij{i,j = 1,2,3) 

^In asymptotically non-free gauge theories, it is easier to have strong convergence, since the gauge 
couplings become large at high energy and b is negative [22]. 
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with the same properties, and extend the Yukawa interactions as 



W ^ YijQiUjHij. 

ij=l,2,3 



(24) 



In the SU(5) GUT, we may take Q = u = 10 and H=H(5). In the supersymmetric 
theories the beta functions for the Yukawa couphngs Yij are written down in terms of the 
anomalous dimensions of Qi, Ui and Hij. This is due to the so-called non-renormalization 
of superpotentials. Explicitly these anomalous dimensions may be written down as 



cqag 



(25) 
(26) 
(27) 



where we defined ay-j = lYij]"^ / (Stt^) as well. Note that the above interactions in Eq. (24) 



induce anomalous dimensions only in the diagonal elements, which are given above. Then 
the RG equation for ay-j is simply given by 



a 



yij- 



(28) 



It is easy to see existence of a non-trivial fixed point by using the beta function given 
by Eq. (28). The RG equations for Xij = CYy^j/cXg are written down as 



d/i 



(6 - c) + {yq. + %. + 7h,, 



where we have defined c — cq + Cu + ch and also 

iQi = aQ{xii + Xi2 + Xi^), 

%i = a^ixij + X2j + X3j), 



(29) 



(30) 
(31) 
(32) 



Therefore the condition for the non-trivial fixed point is that the combinations in the 
brace of the Eq. (29) vanish for all sets of (i, j). The solution is unique and is found to be 



X. 



X 



3a 



(33) 



where a = ag + a„ + an- 

The non-trivial but important matter is whether this fixed point is really IR attractive 
or not. This may be seen by linear analysis of the RG equations around the fixed point 
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(33). The deviations from the fixed point Axjj satisfy the following equation; 



dAxij 
d/i 
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ag 


a' 


ag 




















au 


ag 


ag 


a' 


) 



AXij , 



(34) 



where a' = ag + + 3a//- The eigenvalues of this matrix are found to be 3a//, 3a//, 3a//, 
3a//, 3(aQ + a//), 3(aQ + a//), 3(a„ + a//), 3(a„ + a//), 3a. It should be noted that all of 
these are positive, which ensures IR attractiveness of the fixed point. Thus all Yukawa 
couplings may be aligned to the same value by renormalization effect. ^ 

It is possible also to reduce the number of Higgs upto three by assuming a symmetry 
as follows. We assign the Zs-charge qi = 2n/3i to {Qi,Ui, Hi){i = 1,2,3) respectively. 
Then the symmetry restricts the Yukawa interactions to 



i,i=l,2,3 



(35) 



where the indexes are defined modulo 3. With these interactions the anomalous dimen- 
sions are given to be diagonal. Explicitly they are found to be 



7Qi = 

'■1! = 

The RG equations for 



au{cx.yii + (^y2i + ^V^i) ~ Cu<^g , 
aH{,ayH2-i) + "2/2(l-i) + "2/3(3-i)) ~ ^HOtg 

^vij/^g given by 



1^- 



djj, 



{b-c) + (7Q. + % + 7ff(3_,_,)) 



(36) 
(37) 

(38) 



(39) 



where c 



cg + Cu + Ch and 7 = 7/03 again. It is immediately seen that couplings 
(c — b)/3a (a = ag + a„ + a//) give a set of fixed point solutions. However it is 
found that this is not a unique solution. The linear perturbation around the fixed point 



X, 



X 



^To be precise, what is aligned at low energy is not a Yukawa coupling but it's absolute value. The 
complex phase is not controlled by the dynamics. In this article, we treat the Yukawa couplings as if real, 
and do not discuss the complex phases. 
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satisfies the differential equation, 



dAxij 
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J 



(40) 



where a = aq + + an- The eigenvalues of this matrix are found to be 0, 0, 2aQ, 2aQ, 
2a„, 2au-i 2ajf, 2a/f, 2a. Existence of two zero mode, namely constant modes, means that 
the fixed point solution found above is not IR attractive. In other words, the Yukawa 
couplings obtained at low energy are dependent on their initial values, which are now 
supposed to be disordered. 

However further flavor symmetry may relieve us of this problem. Let us assume e.g. 
a discrete symmetry among matters. 



Ui 



Ui 



Hi Hi. 



(41) 



Then this symmetry imposes additional relations for the Yukawa couplings, which are 
given by 

(42) 



Then it is found that the zero modes in Eq. (40) are just forbidden by this discrete 
symmetry. Also linear perturbation shows us that the fixed point is restricted to a unique 
one and is IR attractive. Thus we have seen that the number of Higgs fields may be 
reduced by assuming some flavor structures in the fundamental theory. However we do 
not consider this possibility more in this article. 



3.3 Flavor democratic Higgs 

So far we have seen that the gauge dynamics may align the Yukawa couplings to the 
same value at IR, once many Higgs flelds are introduced. This does not lead us to a 
democratic mass matrix immediately. Only one mode composed of many Higgs flelds, 
which is identified with the Higgs field in the MSSM, should be massless and all others 
should be superheavy and decoupled at low energy. Moreover every fundamental Higgs 
field must contain the massless mode by the same amount, otherwise the Yukawa couplings 
are deformed from the democratic ones. 

In practice, we may construct the mass terms for the Higgs fields which ensure the 
above properties [18]. We return to the 9 Higgs model given by Eq. (24). In the case of 
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SU(5) GUT, an example of the superpotential for the Higgs fields may be given by 
W = M^(i/(5)(,+i), - i/(5),,)(i/(5)(,+i), - H{5),,) 

i 

+M^(//(5),(,-+i) - //(5),,)(//(5),(,-+i) - H{5),j). (43) 
j 

Then the mass terms for the Higgs scalar fields, 

Y: |i^(5)(.+i), - i/(5),,f + M'Y1 1^(5)^0+1) - H{5),,\' + ■■■, (44) 

i j 

give rise to just one massless mode. Explicitly the massless mode H is given by 

H^WHij. (45) 

Presence of the massless mode is ensured by the shift symmetry with transformations ^, 

H{5),,^C, H{5),,^C, (46) 

where C and C are constants. Other modes acquire mass of order M, which gives decou- 
pling scale. In the later section, we shall take M to be the GUT scale. 

This massless modes composed of H{5) and H(5) arc identified with Higgs fields if" 
and H*^ in the MSSM respectively. It is important to note that the Higgs field H^j in the 
both sector contains the zero mode H with same factor 

H,j^^H + ..., (47) 

for all Therefore the Yukawa matrix appearing in the MSSM is given to be demo- 

cratic one, once the Yukawa couplings have been aligned enough by the strong gauge 
dynamics. 

It may be said that the basic policy in our consideration is flavor democracy. We 
do not assume any hierarchical structures in the initial Yukawa couplings. In addition 
the superpotential for Higgs field does not have fiavor difference. Therefore, there is no 
essential difference among flavors apriori. 

4 A model based on the SU{b) x SU{b) GUT 

Next we consider the mass hierarchy between the flrst two generations and the third 
generation. The small masses for the first and the second generations are generated by 
slight deviations from the rigid democratic mass matrix. From our standpoint of fiavor 
democracy, these deviations should be traces of disorder in the initial couplings. Mass 

^This structure is common to the deconstruction models[23], where the massless mode is generated as 
a Nambu-Goldstone boson. 
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ratio between up and top quarks is more than 1 : 10^. On the other hand, ratio between 
the GUT scale and the Planck scale, which we suppose to be the fundamental scale, is 
only about 1 : 10^. Therefore we need very strong gauge interaction in order to achieve 
such a hierarchy. We may estimate roughly strength of the gauge coupling by assuming 
that it is not running. It is seen from Eq. (21) that the gauge coupling must satisfies 

(47r)2 ^ 6{c-b) 

This condition implies that the gauge coupling need to be non-perturbatively large. 

Naively the gauge coupling constant of the SU(5) GUT is supposed to be rather 
weak, since unification of the three gauge couplings in the MSSM occurs at rather weak 
coupling regime. Therefore it seems unlikely to make the gauge coupling of GUT so 
strong. Actually this is not the case. The gauge couphng unification is not destroyed in 
the presence of any extra heavy fields belonging to SU(5) multiplets. However corrections 
by the extra matters enlarge the unified gauge coupling. Thus one way to realize strong 
unification is to assume a suitable number of extra heavy fields in the MSSM [18] . 

Here we consider another scenario which may be more attractive ^. Suppose that the 
gauge group is given by a product of two SU{5) groups and is broken spontaneously to 
their diagonal subgroup at the unification scale. We represent the gauge couplings for 
the two groups SU{5y and SU{5)" by g' and g" respectively. Here we assume that all of 
the matter fields as well as the Higgs fields are charged under SU{5y but singlet under 
SU (5)" and, moreover, that that g' is non-perturbatively strong but g" is weak. Then RG 
for the Yukawa couplings is driven by the SU{5)' gauge interaction, and may be attracted 
to their IR fixed point very rapidly. Also the gauge coupling g for the diagonal subgroup, 
which is given hj 1/g'^ = 1/g''^ + given weak. 

The superpotential to be considered here is the same as that given in Eq. (16) except 
for the Higgs fields extended to multi ones; 

W = 10a0jH{5yj + Yi 10,5,i7(5),,- + Y^^ 5,l,i7(5),,-. (49) 

To be explicit, the anomalous dimensions for the matter fields and the Higgs fields are 
given at one-loop level as follows; 

35 3 3 

710. = -^«; + 3E< + 2E<' (50) 

^ j=l 3=1 

94 3 3 

75, = -^«; + 4E<i + E<i' (51) 

lu = 5E<„ (52) 

^One may expect also that power law running of gauge coupling in extra dimensions offers us models 
with strong convergence as well. However we do not pursuit for this direction due to a problem discussed 
in section 7. 
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(48) 



24 

7H(5),, = -y"/ + 6ay"^- + (53) 
24 

lH(p),, = -ya; + 4<., (54) 

where ag = and a^"^. = |Fj^p/87r^ for a = u,d, u. The RG equations for the Yukawa 

couphngs Y^^Yf, and Y^, are written down in terms of these anomalous dimensions as 



^^liT " (710, +710, +7f/(5)i,Ja4, (55) 
= (710, +75, +7^/(5),,) "4, (^^) 

^~d^ = (75,+7i, +7/f(5),,ja4- (57) 

In the next section we consider a dynamical mechanism compatible with large mixing 
angles in the lepton sector. There all of the neutrino Yukawa couplings are made 
decrease to very small values. With this reason, we shall consider RG behavior of other 
Yukawa couplings by neglecting the neutrino Yukawa couplings in this section. 

Now two kinds of Yukawa couplings ay^ and Uy'^ are coupled with each other in the 
RG equations. Regardless of this complication, it is found that a non-trivial fixed point 
exists and is IR attractive. The coupling ratios a;^- = O-y^Joig and xf^ = oiy'ljoig satisfy 
the following equations, 

= -- + 6 + ^(3xr, + 3xJ, + 24 + 34)+64U;4, (58) 
A* \ ^ fe=i / 

^ L^ + ^, + ^(3^« +24 + 44)+44)a;4. (59) 

Then it is straightforward to find the non-trivial fixed point solution, which turns out to 
be 

4 = = ^^T5^ ~ 0-^3 - 0.026, (60) 

4 = x^* = ^^^=^ ~ 0.55 - 0.046. (61) 

Thus the Yukawa couphngs are fixed at the GUT scale, once h is given. In any case, the 
Yukawa couplings become non-perturbatively large accompanied with the gauge coupling. 
It is a rather tedious problem to verify the IR attractive nature of this fixed point by linear 
perturbation. However it is quite obvious to see it by solving the differential equations 
given by (58) and (59) numerically. In Fig. 1, the fiow lines for x^- and xfj are shown 
in the case of 6 = 0, aj ~ 1.0. Their initial conditions are chosen at random just for 
demonstration. It is seen that both couplings converge their fixed point values very 
rapidly. 
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Figure 1: RG running of x^- and xfj for a sample of the initial couplings chosen at random. 
We assumed also 6 = 0, ag = 1.0. The parameter t represent the renormalization scale fi 
by using ln(/x/Mpi), where Mpi is the Plank scale. 



People have studied GUT models based on a product group for the purpose of avoiding 
the doublet-triplet splitting problem [19, 20]. Recently models based on G = SU{5)' x 
SU{5)" also have been proposed in this context [20], where gauge charge assignment for 
Higgs fields is different from the above setting. It is assumed that H{5) carries (1,5) 
charges for the product group SU{5y x SU{5)", while H{5) carries (5,1) as before. It 
is noted that the down-type Yukawa interactions in the superpotential given by Eq. (49) 
are not invariant any more. In order to generate the down-type Yukawa interactions, 
we introduce another field S belonging to (5, 5). We suppose also that the spontaneous 
symmetry breaking S'f/(5)' x SU{5)" 5'f/(5)diag is induced by a vacuum expectation 
value of S. Then an invariant, but non-renormalizable, term 

^ 10,5,if(5)i,S (62) 

may generate the down-type Yukawa couplings effectively after symmetry breaking. It 
seems rather difficult to discuss renormalization involved with non-renormalizable opera- 
tors like this. However we may also consider the corrections at the broken vacuum, where 
these operators plays a role of the Yukawa interactions instead. Therefore it is expected 
that the above analysis for the IR fixed point is equally applied to this model. 

As a general property of the present mechanism, the Yukawa couplings for the third 
generations are given to be rather large at the GUT scale. The masses for top, bottom 
and tau are given in terms of the Yukawa couplings obtained after diagonalization. It is 
easy to find these Yukawa couplings renormalized at O(100)GeV by solving RG equations 
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for the MSSM. The neutrino Yukawa couphngs may well be neglected in this analysis. 
Then the low energy couplings are determined with respect to the strong gauge coupling 
g' and the model parameter h. However it is found that the resultant Yukawa couplings 
are almost insensitive to these parameters. This is because the initial Yukawa couplings 
are fairly large and such flows converge to the so-called quasi-fixed point [21]. To be 
exphcit, the low energy couplings are found to be Yt — 1.00, = 0.94 and y,- = 0.62 in 
the case of ag = q'^/'&'k'^ — 1.0 and 6 = 0. It is seen first that tan/? should be large as 
0(50), to explain the mass ratio mt/m^. since Yt and Yf, appear to be almost the same. 
Therefore mass of top quark is predicted to be the same as the vacuum expectation value 
of the neutral Higgs 174GeV, which shows a very good agreement with observation. We 
may also predict mass ratio of bottom and tau, m(,/m^. The above couplings bring us to 
mij/mj- — 1.52, which is slightly smaller than the expected value.^ 

5 Neutrino sector 

In the previous section we have seen that the Yukawa couplings arc aligned to the demo- 
cratic forms very well. The deviations from the rigid democratic form are responsible for 
masses of the first two generations and mixing angles among quarks. However the mixing 
matrix becomes necessarily close to an identity matrix for any small deviations. If we take 
the neutrino Yukawa couplings, ignored in the previous section, into account, then these 
couplings are found to be attracted to an IR fixed point by strong dynamics as well as 
other Yukawa couplings. Obviously this is incompatible with the observed large mixing 
angles for leptons. 

In this respect, it is phenomenologically favorable for the neutrino Yukawa matrix to 
be nearly diagonal at low energy. Suppose that the neutrino Yukawa couplings do not 
have a non-trivial fixed point with some reasons. Then these couplings are not aligned to 
the democratic form and become dependent on their initial couplings at the Planck scale. 
Also we may find a sizable parameter region for the initial couplings consistent with the 
observed neutrino mass and mixings. 

In practice such a situation turns out to be realized with a simple assumption, since the 
right-handed neutrino is a singlet of the 5'f/(5). In general, GUT models may contain a 
pair of superheavy vector-like fields and so on. Then the right-handed neutrino is allowed 
to have additional Yukawa interactions with them These Yukawa couplings are driven 
to be very large like others by the strong gauge interactions. Consequently anomalous 
dimension of the right-handed neutrino is enhanced significantly. We note that the IR 
fixed point stands upon balance of negative contribution by gauge interactions and positive 
contribution by Yukawa interactions. Presence of extra Yukawa couplings may destroy 
this balance and take the IR fixed point away from the neutrino Yukawa couplings in the 
end. Explicitly the coefficient c in the fixed point equation (29) is made negative by the 

^In the case with large tan /3, we have to take into account large SUSY threshold corrections on the 
bottom mass [24]. 

^°R-parities for the extra fields may be assigned properly. 
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additional radiative corrections. Thus the remarkable difference in the mixing matrices 
of quark and lepton sectors may attribute to the dynamics of the right-handed neutrino. 

We shall demonstrate the above mechanism by examining the RG equations only for 
the neutrino Yukawa couplings. As a toy model with a single flavor, let us consider the 
super potential given by 

51 H{5) + k1 $ l>, (63) 

where the extra matter $ and i» belong to SU (5) representations R and R respectively. 
The RG equations for ay — |yp/87r^, and q;^ = |«;p/87r^ are found to be 

^~d^ = [lay + Ra^-AC2{5)ag]ay, (64) 

dtoi 

= [hay + {R + 2)a^-AC2{K)ag]a^, (65) 

where C2 denotes the Gasimir index of each representation. Also these are rewritten into 
equations in terms of Xy — ay jag and x^ — a^/agas 

dx 

H-^ = [7xy + Rx^-4C2i5) + b]agXy, (66) 
dx 

= [bXy + {R+2)x^-AC2{B) + h]agX^. (67) 

There are four flxed points in the coupling space of {xy, x,^), three of which are immediately 
seen from the above equations as 

{xy,x,) = (0,0), (68) 

■^^^'^'-^ol. (69) 

The first two fixed points are what we have been discussing so far and the third one is a 
new entry. It is noted that the fourth one appears in the region of negative ay for a large 
representation R Here we pay attention to such a case. 

The RG flow diagram for {xy, x,^) is shown in Fig. 2 in the case of = and R = adj 
as an example. The marked points A, B and G stands for the flxed points given by (68), 
(69) and (70) respectively. It is seen that the fixed point B is now unstable towards 
direction of the new coupling, while C turns out to be IR attractive instead. Thus the 
neutrino Yukawa coupling is found to decrease at low energy in the presence of k. 

In Fig. 3 evolution of neutrino Yukawa couphng with respect to the scale is shown. 
Each line corresponds with a RG flow presented in Fig. 2. It is found that the Yukawa 




^^This point may not be unphysical with taking complex phases of the Yukawa coupling into consider- 
ations. However this fixed point is not IR attractive anyway. 
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A B 

Figure 2: RG flows for {xy, x^) in tlie case of 6 = and R = adj. Tfie arrows show the 
direction to lower energy scale. The circles at A, B and C stand for the fixed points given 
by (68), (68) and (70) respectively. 



couplings decrease monotonically at low energy and are not aligned. It is expected that 
the neutrino Yukawa couplings for multi-flavors also show the same aspect of these flows, 
even though other Yukawa couplings Y"^ and are aligned to the fixed point. We now 
add the following terms, 

^5iljH{5)ij + Hih<!> $ + Mnijhlj, (71) 

to the superpotential given by Eq. (49). It is easily found that ratio of the neutrino 
Yukawa couplings follows the RG equation given by 

where the Yukawa couplings other than ay'^ are assumed to be their fixed point values. 
Since the neutrino Yukawa couplings are suppressed as seen above, the right-hand side of 
(72) becomes very small quickly. This shows that ratio of the neutrino couplings is almost 
unchanged at low energy region. 

In the Frogatt-Nielsen mechanism [6], the neutrino Yukawa couplings are constrained 
upto 0(1) coefficient, and, therefore, the masses and mixing angles for neutrinos are not 
predicted. However we may find a sizable probabihty for the couplings distributed in 
anarchy to be compatible with the observed masses and mixing angles [25]. 
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Figure 3: Running of Xy corresponding with various flow lines shown in Fig. 2 with respect 
to the scale fJ. (t = ln(/x/Mpi)). 



Now if three of nine initial neutrino Yukawa couplings happen to be larger than others 
to some extent, then they may be dominant in the low energy matrix also. Since there 
is no specific flavor basis, we may regard these couplings as the diagonal elements. Thus 
it is not so special that the neutrino Yukawa coupling matrix appears close to diagonal, 
and therefore the large lepton mixing angles are described well. We leave explicit survey 
of such a parameter space for future works. 

In our mechanism the neutrino Yukawa couplings are largely suppressed as seen in 
Fig. 3. However this does not imply that the see-saw mechanism offers us very tiny 
neutrino masses. Note that this suppression occurs because of enhancement of anomalous 
dimensions of the right-handed neutrinos. Then their Majorana masses Mn are also 
suppressed by the anomalous dimensions. Consequently it is found that the neutrino 
masses themselves are not suppressed through this mechanism. It is noted that the 
masses for the right-handed neutrinos are expected to appear at the intermediate scale in 
our scenario, even though the bare parameter for Mr is set to be the GUT scale. 

6 Setup for the realistic mass matrices 

How can we describe mass differences among the first and the second generations at all? 
In the phenomenological approach, the two flavor symmetry breaking parameters e and S 
are taken to be hierarchical as given in Eq. (6). Besides it is necessary for these breaking 
parameters to be diagonal elements, which is not explained from the standpoint of flavor 
symmetry either. 

Now the problem to be considered is how to realize such a specific form of deviations 
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from the democratic couplings in our dynamical framework. The first possibility is to 
generate such deviations away from the fixed point by assuming proper initial couplings. 
We consider this by examining the RG equations given by Eq. (29) again. Suppose, e.g. 
that one of the Yukawa couplings, say 133, happens to be fairly smaller than the others 
at the Planck scale. Then Xij for i,j — 1,2 converge into the fixed point value rapidly, 
since the beta functions for these couplings do not contain X33. The couplings and 
for i — 1,2 may be affected by ,T33. So we may examine the RG fiows with neglecting 
difi^erences among Xij for i, j = 1,2 in the first step of approximation. In this situation 
the deviations from the fixed point Axij found to satisfy 

[Axi3 - Aa;23] = agX* {a' - ttu) [Axis - Ax23] (73) 

/x^ [Axsi - AX32] = agX*{a' - ag) [Ax^i - AX32] ■ (74) 

These equations tell us that these differences shrink rapidly with scaling down and even- 
tually the couplings satisfy the relations 

3:^13 = X23, X31 = X32- (75) 

Moreover differences among Xij for i,j = 1,2, which is ignored in the above argument, 
are found to satisfy 

H-^ [Axu - Ax2i] = fjL-^ [Axi2 - AX22] = agX*aQ [Axi3 - AX23] , (76) 

fjL-^ [Axu - AX12] = //^ [Ax2i - AX22] = agX*au [AX31 - AX32] . (77) 

Therefore it is seen that these differences are not affected from large deviation of 2:33 and 
are kept small. In Fig. 4 a sample of RG flows for the couplings Xij obtained by numerical 
analysis are demonstrated. The bold line shows flow of X33, whose initial coupling is 
assumed to be relatively smaller than others. Note that Yukawa coupling Y33 itself is not 
so apart from others. The dashed lines represent flow of Xi^ and Xsi for i = 1,2, while 
the dotted lines represent flow of Xij for i,j = 1,2. Indeed behavior of these RG flows 
supports the above argument. 

Similar forms of Yukawa coupling matrices may be obtained by introducing a weak 
interaction breaking the flavor democracy. For example, we assume only one of the Higgs 
flelds, say ^^33, has an extra superpotential Ai/33<I>i<I>2 with very small coupling A. We 
consider the linear perturbation around the fixed point. The extra Yukawa interaction 
affects only on the beta function of Aa;33. To be explicit, the RG equations given in 
Eq. (34) are modified to 

r/Ar ■ 

^^if ^ ii^^^h + '^A<^i3<^,-3] , (78) 
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Xij 

1 r 




-1.5 -1.25 -1 -0.75 -0.5 -0.25 

Figure 4: Convergence behavior of Xij to the IR fixed point in the case of a proper initial 
couphngs. The parameter t denotes ln(/i/Mpi). The bold line shows fiow of X33, whose 
initial coupling is assumed to be rather smaller than others. The dashed lines represent 
fiow of Xj3 and xsi for i = 1,2, while the dotted lines represent fiow of Xij for i,j = 1, 2. 



where Ai denotes the matrix given in Eq. (34) and 6\agX* ~ |Ap/87r^. When differences 
among the initial Yukawa couplings happen to be rather small already, the same equations 
given by Eqs. (73), (74), (76) and (77) hold as well. 

In the end, the Yukawa couplings turn out to be represented by the form of 



Y = yo 



( l + 0(e) l + 0(e) 1 + 5' 

1 + 0(e) 1 + 0(e) 1 + 5' I , (79) 
V 1 + 5" 1 + 5" 1 + 5 



at low energy scale. Here e stands for the incompleteness of strong convergence to the 
fixed point, which can be as small as 10~^. Other parameters representing deviations from 
the fixed point 5, 5' and 5" are supposed to be of the same order. Though these values 
are dependent on the initial Yukawa couplings, we may well assume e <C 5 <C 1. After 
transforming this matrix Y by using the diagonalizing matrix given by Eq. (2), we obtain 

/0(e) 0(e) 0{5) \ 
A'YA = yA 0(e) 0{5) 0{5) . (80) 
V 0{5) 0{5) 3 + 0(5) ; 

This is very similar to the Fritzsch type matrix and, therefore, offers us a phenomenologi- 
cally favorable pattern for realistic masses and mixing angles. The ratio mass eigenvalues 
is given by 0(e) : 0(5) : 0(1). We do not study which initial couplings really give the 
viable mass matrices explicitly, since there are large model dependence. 
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7 Conclusions and discussions 



The approach of democratic mass matrix seems to be attractive enough phenomenolog- 
ically. The notion of flavor democracy offers us an interesting viewpoint distinct from 
other approach. Above all it is noted that the lepton mixing angles, bi-large mixings 
and small 1/^3, are well described by assuming a democratic mass matrix for the charged 
leptons. 

However it is also true that these phenomenologically favorable matrices are not ex- 
plained by flavor symmetries. In this paper we considered the models in which the demo- 
cratic type of Yukawa couphngs are reahzed as an IR attractive fixed point [18]. For this 
purpose we introduced multi Higgs fields, from which only one massless mode is allowed 
by a special type of mass terms. The large mass hierarchies are achieved easily by using 
strong gauge interaction for quarks and leptons as well as the Higgs fields. To be explicit, 
we considered an SU (5) x SU (5) GUT model, in which one of the gauge couplings is 
non-perturbatively strong. 

Indeed the setup of the Higgs fields may be somewhat artificial. However we found 
some benefits in this scenario. The democratic Yukawa couplings are realized even in 
GUT models without any adhoc assumptions. Hence small mixing angles as well as large 
mass hierarchies in the quark sector are naturally obtained. Especially the large lepton 
mixing angles can be explained through extra interactions of the right-handed neutrinos. 
It is seen that simple setting of the initial couplings may lead to low energy Yukawa 
couplings of the Fritzsch type. Also this kind of scenarios predict top quark mass in the 
good agreement with experiment and large tan f3. 

In this article we considered only supersymmetric models. However the PR fixed point 
exists also in non-supersymmetric theories as well. Therefore the democratic mass matri- 
ces may be realized also in some non-supersymmetric models, though the beta functions 
for Yukawa couplings are complicated. 

In section 3 we mentioned that the PR fixed point appears also in extra dimensions 
and shows rather strong convergence in general. However naive extension with gauge and 
Higgs in the bulk of the extra dimensions does not work. This is because anomalous 
dimensions of the Higgs fields vanish due to effective N=2 supersymmetry, and , therefore, 
the democratic type of fixed point is not realized. This would not deny any possibility 
for extra dimensional models, and we leave it for future study. 

Lastly some comments on soft supersymmetry breaking parameters would be in order. 
It has been known [27, 28] that the PR fixed point induces pecuhar relations among soft 
parameters also. Above all A-parameters are found to be aligned in the present scenarios 
with the same dynamics for the Yukawa coupling alignment. The soft scalar masses for 
quarks and leptons become fiavor universal at low energy due to large corrections by the 
democratic Yukawa couplings as well as the gauge couplings. These properties are very 
desirable for the flavor problem in supersymmetric extensions. 

However the strong gauge dynamics enhances the soft scalar masses and also the 

^^The various flavor violating processes have been examined in the phenomenological approach of the 
democratic mass matrices [26]. 
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A-parameters to be comparable with the strongly coupled gaugino mass. Therefore the 
gaugino mass in the strongly coupled sector must be much smaller than the MSSM gaugino 
masses. In this respect, extension to extra dimensions will be considerable, since FCNC 
and also CP problems may be solved or ameliorated without assuming such a special case 
[29] . Otherwise it would be more natural for our scenarios that supersymmetry breakings 
are mediated below the GUT scale e.g. by gauge mediation mechanism. Further studies 
of the supersymmetry breaking parameters in our models will be discussed elsewhere. 

Acknowledgements 

T. K. is supported in part by the Grants- in- Aid for Scientific Research (No. 16028211) and 
the Grant-in- Aid for the 21st Century COE "The Center for Diversity and Universality 
in Physics" from the Ministry of Education, Science, Sports and Culture, Japan. H. T. 
are supported in part by the Grants-in-Aid for Scientific Research (No. 13135210) from 
the Ministry of Education, Science, Sports and Culture, Japan. 

References 

[1] Super-Kamiokande Collaboration, Y. Fukuda et al, Phys. Lett. B 467, 185 (1999); 
S. Fukuda et al, Phys. Rev. Lett. 85, 3999 (2000); Phys. Rev. Lett. 86, 5651 (2001); 
Phys. Rev. Lett. 86, 5656 (2001). 

[2] SNO Collaboration, Q. R. Ahmad et al, Phys. Rev. Lett. 87, 071301 (2001); 89, 
011301 (2002); 89, 011302 (2002). 

[3] KamLAND Collaboration, K. Eguchi et al, Phys. Rev. Lett. 90, 021802 (2003). 

[4] K2K Collaboration, M. H. Ahn et al, Phys. Rev. Lett. 90, 041801 (2003). 

[5] For recent reviews, see for example: H. Fritzsch and Z. Z. Xing, 
Prog. Part. Nucl. Phys. 45, 1 (2000); S. F. King, Rept. Prog. Phys. 67, 107 (2004); 
G. AltareUi and F. Ferugho, hep-ph/0405048. 

[6] C. D. Froggatt and H. B. Nielsen, Nucl. Phys. B 147, 277 (1979). 

[7] See e.g., S. F. King and G. G. Ross, Phys. Lett. B 520, 243 (2001); Phys. Lett. B 
574, 239 (2003). 

[8] H. Fritzsch, Nucl. Phys. B 155, 189 (1979); Y. Koide, Phys. Rev. D 28, 252 (1983); 
Phys. Rev. D 39, 1391 (1989). 

[9] H. Fritzsch and Z. Z. Xing, Phys. Lett. B 372, 265 (1996); Phys. Lett. B 440, 313 
(1998); Phys. Rev. D 61, 073016 (2000); Phys. Lett. B 598, 237 (2004). 



23 



[10] M. Tanimoto, Phys. Lett. B 483, 417 (2000); M. Fukugita, M. Tanimoto and 
T. Yanagida, Phys. Rev. D 57, 4429 (1998). 

[11] M. Fukugita, M. Tanimoto and T. Yanagida, Phys. Rev. D 57, 4429 (1998). 

[12] M. Tanimoto, T. Watari, T. Yanagida, Phys. Lett. B461 (1999) 345; E. Kh. Akhme- 
dov, G. C. Branco, F. R. Joaquim and J. I. Silva-Marcos, Phys. Lett. B498 (2001) 
237; G. C. Branco and J. L Silva-Marcos, Phys. Lett. B526 (2002) 104; T. Watari 
and T. Yanagida, Phys. Lett. B544 (2002) 167; Q. Shafi and Z. Tavartkiladze, hep- 
ph/0401235. 

[13] G. C. Branco and J. I. Silva-Marcos, Phys. Lett. B 359, 166 (1995); G. C. Branco, 
M. N. Rebelo and J. I. Silva-Marcos, Phys. Lett. B 428, 136 (1998); Phys. Rev. D 62, 
073004 (2000); E. K. Akhmedov, G. C. Branco, F. R. Joaquim and J. I. Silva-Marcos, 
Phys. Lett. B 498, 237 (2001). 

[14] A. E. Nelson and M. J. Strassler, JHEP 0009, 030 (2000); JHEP 0207, 021 
(2002); T. Kobayashi and H. Tcrao, Phys. Rev. D 64, 075003 (2001); T. Kobayashi, 
H. Nakano and H. Terao Phys. Rev. D 65, 015006 (2002); T. Kobayashi, H. Nakano, 
T. Noguchi and H. Terao Phys. Rev. D 66, 095011 (2002). 

[15] M. Bando, T. Kobayashi, T. Noguchi, K. Yoshioka, Phys. Lett. B 480, 187 (2000); 
Phys. Rev. D 63, 113017 (2001). 

[16] B. Pendelton and G. G. Ross, Phys. Lett. B 98 291 (1981). 

[17] G. G. Ross, Phys. Lett. B 364, 216 (1995); B. C. AUanach and S. F. King, 
Phys. Lett. B 407, 124 (1997). 

[18] S. A. Abel and S. F. King, Phys. Lett. B 435, 73 (1998). 

[19] R. Barbieri, G. Dvah and A. Strumia, Phys. Lett. B 333, 79 (1994); T. Yanagida, 
Phys. Lett. B 344, 211 (1995); T. Hotta, Izawa K.-L, T. Yanagida, Phys. Rev. D 53, 
3913 (1996). 

[20] E. Witten, hep-ph/0201018; M. Dine, Y. Nir and Y. Shadmi, Phys. Rev. D 66 115001 
(2002). 

[21] C. T. Hill, Phys. Rev. D 24, 691 (1981). 

[22] M. Bando, J. Sato and K. Yoshioka, Prog. Theor. Phys. 98, 169 (1997); 
Prog. Theor. Phys. 100, 797 (1998). 

[23] N. Arkani-Hamed, A. G. Cohen and H. Georgi, Phys. Rev. Lett. 86, 4757 (2001); 
Phys. Lett. B 513, 232 (2001); C. T. Hill, S. Pokorski and J. Wang, Phys. Rev. D 
64, 105005 (2001). 



24 



[24] L. J. Hall, R. Rattazzi and U. Sarid, Phys. Rev. D 50, 7048 (1994); R. Hempfling, 
Phys. Rev. D 49, 6168 (1994). 

[25] L. Hall, H. Murayama and N. Weiner, Phys. Rev. Lett. 84, 2572 (2000); N. Haba and 
H. Murayama, Phys. Rev. D 63, 053010 (2001); A. de Gouvea and H. Murayama, 
Phys. Lett. B 573, 94 (2003). 

[26] K. Hamaguchi, M. Kakizaki and M. Yamaguchi, Phys. Rev. D 68 056007 (2003). 

[27] M. Lanzagorta and G. G. Ross, Phys. Lett. B 364, 163 (1995); S. F. King and 

G. G. Ross, Nucl. Phys. B 530, 3 (1998). 

[28] T. Kobayashi and K. Yoshioka, Phys. Rev. D 62, 115003 (2000); T. Kobayashi and 

H. Terao, Phys. Lett. B 489, 233 (2000). 

[29] J. Kubo and H. Terao, Phys. Rev. D 66, 116003 (2002); Y. Kajiyama, J. Kubo and 
H. Terao, Phys. Rev. D 69, 116006 (2004); K-Y. Choi, Y. Kajiyama, H. M. Lee and 
J. Kubo, Phys. Rev. D 70, 055004 (2004). 



25 



